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Abstract 

Exploiting the cone structure of the set of unnormahzed mixed quantum states, we offer an 
approach to detect separabihty independently of the dimensions of the subsystems. We show that 
any mixed quantum state can be decomposed as p = (1 — A)Cp + XEp, where Cp is a separable 
matrix whose rank equals that of p and the rank of Ep is strictly lower than that of p. With 
the simple choice Cp = Mi M2 we have a necessary condition of separability in terms of A, 
which is also sufficient if the rank of Ep equals 1. We give a first extension of this result to detect 
genuine entanglement in multipartite states and show a natural connection between the multipartite 
separability problem and the classification of pure states under stochastic local operations and 
classical communication (SLOCC). We argue that this approach is not exhausted with the first 
simple choices included herein. 
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I. INTRODUCTION 



"Entanglement is [. . . ] the characteristic trait of quantum mechanics" [l| . This state- 
ment remains vahd, even stronger, more than 70 years after its formulation. Nowadays 
entanglement is believed to lie behind most quantum phenomena and is still considered the 
source of conceptual challenges ^|. In a few words, entanglement stands for the presence of 
quantum correlations, i.e. statistical correlations which cannot be reproduced using classical 
probability theory. Despite its importance, we are still far from a full comprehension of 
this fundamental feature of Nature, not to mention the fact that the problem of discern- 
ing whether a given composite quantum system portrays entanglement has only a partial 
solution. By and large, a complete understanding will follow only after apprehending the 
notion of quantum correlation both qualitatively and quantitatively for all sort of compos- 
ite systems, either bipartite or multipartite and either with infinite or with finite dimensions. 



Today a big gap exists between our understanding of bipartite and multipartite entangle- 
ment, probably as big as radical the difference is between having two subsystems and having 
more than two. This gap is also present between finite- and infinite-dimensional systems. 
These latter have been analyzed mainly in the Gaussian domain (see e.g. jsj). In the 
following we will concentrate upon finite-dimensional systems. From a broad mathematical 
point of view two approaches can be followed to discern whether a given quantum state is 
entangled or not. On the one hand we can adopt an analytical point of view and try to 
find operational criteria to settle the question upon the quantum state itself. As prominent 
instances we can cite the Schmidt decomposition (only valid for pure states of any dimen- 
sions) y , the PPT criterion (valid as a necessary and sufficient condition for 2 (g) 2 or 2 ® 3 
and a necessary condition for the rest of bipartite systems) 5|-[3], the range criterion (a nec- 
essary condition for all bipartite systems) [8| and the computable cross norm or realignment 



criterion (a necessary condition for all bipartite systems) [q], |io|, to name a few (see ll|-|l4| 



for an overview from different perspectives). On the other hand, we can adopt a numerical 



viewpoint, by which the separability problem for a given state is trans^ 
mization problem and then an adequate algorithm is chosen to solve it 



ated into an opti- 



ISMlTj. Several of 



these approaches can be generalized to multipartite systems. A remarkable position in this 
view is played by entanglement witnesses which accept a double role in this division. 



The numerical techniques can be understood to solve the question, but, in our view, 
this cannot be completely satisfactory inasmuch as a numerical solution to a second degree 
algebraic equation is not as satisfactory as its usual analytical solution. Furthermore we be- 
lieve that such an analytical comprehension of the separability problem may help in related 
questions such as the distillability problem and possibly open the way to understand 
entanglement in a much larger scale, in particular in quantum many-body systems (see e.g. 



2l|). In this line we find it mandatory to pursue an analytical operational approach valid 



for any composite system in any dimensions. 

Admittedly this is too ambitious a goal to cherish a realistic hope to find a global solution 
in a short term and we must restrict this generality. In this work we present the first steps 
aiming at finding an analytical operational solution to the bipartite entanglement problem 
in any dimensions. Our proposal rests on the well-known structure of the set of Hermitian 



positive semidefinite matrices, namely, the convex cone structure PSBtvi^atj |22|, valid for 
any dimensions A^i and N2. Our bottom line tries both to be independent of the dimensions 
and of the bipartite character of the composite system, seeking a plausible generalization to 
multipartite systems. The price to be paid is that we have not covered all possible quantum 
states yet (see below). However we believe our results can be pushed forward to encompass 
the remaining cases and we hope to do this in a near future. 



Our approach identifies the set of separable matrices SEPTVi^Afa as a subcone of PSD 



and exploits the cone properties of this to write any unnormalized 23| quantum state 
Q G PSD7Vi(g)Ar2 as a convex combination Q = {1 — Xq)Cq + XqEq of a separable matrix Cq 
with the same rank as Q and a matrix Eq in the boundary of PSBtvi^atj- There are similar 



decompositions in the literature |2^. l25j but with a different spirit. Our goal is to determine 
those values of Ag G [0, 1] by which Q is separable and those for which it is entangled. As 
a result, we will also have an explicit decomposition of Q in terms of separable matrices. 
Regretfully for the time being we have only partially succeeded for those matrices whose 
range S = R{Q) can be written as a tensor product S = Si ^ S2. We also point out 
how this formalism admits a generalization to the multipartite case through the SLOCC 



entanglement classification of pure states 



12,l26|,l27|. 



The paper is organized as follows: in section |TT] we review those pertinent cone properties 
of the set of Hermitian positive semidefinite matrices. In section UTT] we give the bottom line 
of our approach, which is applied in section IIVI to full-rank quantum states and in section 
|V]to quantum states lying in product faces. In section |VT] plausible arguments to generalize 
these techniques to the multipartite case are given as well as the natural connection with 
the SLOCC classification of pure states. We conclude with final remarks in section I VI II 



II. THE CONE STRUCTURE OF UNNORMALIZED QUANTUM STATES 

A natural mathematical framework for finite-dimensional quantum states in the light 
of the separability problem is the cone of Hermitian positive semidefinite matrices. We 
review very briefiy the notion of a cone and the properties of PSDjv^^jVj, the cone of these 
matrices. This amounts to relaxing the separability problem to find a decomposition of an 
unnormahzed quantum state Q G PSD^v-^^atj as 



M 



Q = ^\i\xi){xi\®\yi){yil Aj > Vz = 1, 



M. 



(1) 



i=l 



Note that the non-negative parameters Aj are not restricted to satisfy YldLi^i = 1, as 
in the original formulation 28|. We will very briefiy include those notions pertinent to 
the forthcoming approach in order to render this work as self-contained as possible. The 
interested reader may consult the references 



22 



29| for more details. 



Definition II. 1. Let V he a real topological vector space. A cone K in V is a closed subset 
ofV which satisfies (i) if a, /3 > and x, y G then a;x + /3y G K and (ii) Kr\ {—K) = 0. 

The algebraic structure of a cone generalizes in some sense that of a vector space. In this 
line of thought, the analogous roles of the concepts of vector subspace and linear basis are 
played by those of faces and extremals, respectively. 

Definition II. 2. A suhcone F d K is a face of the cone K if 

X G -ft', y — X G -fC, and y G -F imply x G -F. 
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Implicit in this definition is the fact that a face is a cone itself. Let dimF = dim(F — F) 
and $(5) denote the least face of K containing the vector subspace S. 

Definition II. 3. A face F is an extreme ray of the cone K if dim F = 1. A vector x is an 
extremal z/ dim $(sj)an{x}) = 1. 

It will be often useful to identify a cone by its extremals, since any vector in the cone 
can be expressed as a so-called conic combination of extremals (except the extremals them- 
selves), i.e. y = J27=i^i^i' where Aj > and Xj are extremals. Note how close is this 
expression to that of a linear combination in terms of elements of a linear basis. In this 
sense it is written K = cone{xj}tgT, where all vectors xj are extremals and the index set T 
may be non-numerable j30|. 



All these concepts apply readily to the set of Hermitian positive semidefinite matrices 
acting upon the complex vector space C^^ (g) C^^. This set is a cone within the real vector 
space 'HN-i®N2 of Hermitian matrices denoted by PSDAr^07V2- In particular, the extremals are 
matrices of the form A|x)(x|, the extreme rays are subsets of the form {A|x)(x|}a>o and the 
faces are parametrized by the different vector subspaces belonging to ® in such a 
way that given a vector subspace iS, then the face F = is the set of all positive matrices 
unnormalized quantum states) Q such that R{Q) C S, where R{Q) denotes the range of Q 



22|. Note that we can write 



PSDAT^gjAT^ = cone{|x)(x|}|^.)gci®C2 (2) 

Indeed the spectral decomposition of any Hermitian positive semidefinite matrix is an 
example of this fact, since Q = Xlil^^^ K\xi){xi\, where Aj > are the eigenvalues of Q and 
\xi) the corresponding eigenvectors. Notice also that ([2]) expresses the fact that this kind of 
decompositions is valid without necessarily having \xi) _L \xj) for all i ^ j. Furthermore the 
faces, being cones themselves, can also be written as 

Fs = cone{|x)(x|}|^)e5. (3) 

Important characteristics of PSDiVi®Af2 to be exploited later are (see j22| and references 
therein) : 



1. The set of non-trivial faces (faces not being the whole cone PSDAr^07V2 itself) constitutes 
the topological boundary of PSD7v-^0Ar2, which we will denote by SPSDatj^at^ |31[ |. 



2. The preceding result is also valid for each face considered as a cone in the relative 
topology. 

3. FSDtv^^at^ and any of its faces are convex. 

4. PSDat-^^atj is connected and its faces are also connected in their relative topology. 

5. Since every face F = Fg is characterized by a vector subspace S we can attach to each 
face an integer, namely r = dimiS; this is indeed the rank of the greatest orthogonal 
projector Ps = = Pj belonging to Fs, i.e. r = rp^ = dimR{Ps). 

As an example of the advantages of this approach, we offer a one-line proof that the cone 



32| [33|: any face F = Fg^ {dimSk = k) of PSDjy has 



of separable matrices has no facet 

dimension k"^ {22!, thus any subcone of PSBtvi^atj has no facet, let alone the cone SEP; 



III. THE GENERAL IDEA 



The proposal included herein rests upon the two following independent results. We will 
denote by int {A) the topological interior of set A. 

Proposition III.l. Let Q E inti^Fs^^. Then there exists another matrix Cq G inti^Fg^^ 
and a matrix Eg G dF^^ such that 

Q = tr{Q) [(1 - Xq)Cq + XqEq] , (4) 
where Xq G (0, 1) and triCq) = tr^Eq) = 1. 

Proof. With no loss of generality we can suppose that tr((5) = 1 (otherwise repeat the 
forthcoming argument for Let Ai denote the affine subspace of unit-trace Hermitian 

positive semidefinite matrices. Since Q G (int (-P^g)) H Ai there always exists an open 
neighborhood of Q contained in (int (-F^g)) CiAi. Take Cq as any of the matrices belonging 
to this neighborhood. Since (Fsq) H Ai is connected it is possible to construct the closed 
segment (1 — jj)Cq + /iQ, with /x G [0, 1], lying inside (int {Fsq)) flAi for all and to extend 
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the segment (fi > 1) until (1 — ii)Cq + G (OFsq) f^Al. Let /i* be the minimal yU > 1 such 
that (l-;u)CQ + /ig e [dFsQ)nAi. Then = (l-;U*)CQ + /i*g e (SF^^) flAi and, given 
the convexity of Fs^ flAi there always exist a real A G (0, 1) such Q = {1 — \)Cq + XEq. □ 



Note that decomposition (jlj) is similar to others pursued in different approaches |24j. |25|. 
The resemblance arises, in our view, from the convex cone structure of the set of quantum 
states; the difference stems from the diverse alternative we have when choosing either Cq 
or Eq. This decomposition is rather easy to compute from a theoretical point of view: find 
the minimal > 1 such that r(i_^)CQ+^Q < tq. 



The second result we need is the following 
Proposition III. 2. The cone of separable matrices SEPn^,^m2 ^ proper suhcone of 

Proof. This elementarily follows from the fact that the set of extremals {|a;)(a;|®|y)(?/|}|^.)gcivi |y)gc^2 
of SEPat^^tvj is strictly contained in the set of extremals {|-2)}|2>ec^i®c^2 of SEPatji^atj. □ 

The general idea of this approach lies on the choice of Cq in the decomposition (jlj) in the 
light of proposition IIII.2I 

Criterion III.l. Given a quantum state p choose Cp G SEPn^^n^ and let A* be the greatest 
Ap such that (1 — A*)Cp + A*£'p G SEPn^c^)]^^. Then p will be separable if and only if Xp < A*. 

Proof. The result follows readily from the convexity of both SEPjVi^ATa and PSD^Vi^ATa and 
from the fact that SEPat-^^atj C PSDat-^^atj. Since both Cp and (1 — A*)Cp + \*pEp are 
separable, then (1 — Xp)Cp + XpEp is separable whenever Ap G [0, A*], and for Ap > A*, by 
construction, (1 — Xp)Cp + Api?p is not separable. □ 

Note that given a quantum state p this separability criterion is reduced to the face Egp, 
thus not involving quantum states belonging to the rest of the cone. Although not indicated 
in the notation, the A* values depend also on Cp and Ep. Notice also that once the values 
A* are determined, the separability check consists of a simple comparison between the value 
Ap in the decomposition (jlj) and the corresponding A* value. In the subsequent sections we 
will concentrate upon the computation of these values. 
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IV. FULL-RANK QUANTUM STATES 



Firstly we will concentrate upon finding the A* values for those quantum states p having 
full rank. Following proposition II 1 1 . 1 1 we have always the freedom to choose either Cp or Ep. 
Our choice is Cp = Mi ® M2, where Mj G int PSBtv^, that is, = A^j for i = 1, 2. The 
resulting Ep in decomposition (jl]) is not determined a priori, so it could be any matrix with 
rank strictly lower that Ni- N2. We will distinguish between the cases in which = 1 and 
those in which te^ > 1. 

We need previously some notation and some preliminary results. Given a pure quan- 
tum state \z) G C^^ (g) C^^, we define the {Ni x A^2)- dimensional matrix of coefficients in 
the computational basis \z) = XliiJa ^*i*2 ^1^2) by C{z) = [cj^ij], where ii is a row column 
and 22 a column index. Then we can apply any generalized singular value decomposition 
(GSVD) to C{z) so that C{z) = F{z)D{z)G''{z), where F{z) is an A''i-dimensional square 
nonsingular matrix, D{z) is a (A^^i x dimensional diagonal matrix with nonnull entries 
c"o(^) > ...crj,_i(z) > and G{z) is an Ai'2-dimensional square nonsingular matrix (see 
appendix |X] for an overview of the GSVD). For ease of notation we do not specify the 
particular GSVD yielding the generalized singular values Cilz), this being clear from the 
context. When the proper SVD is indeed applied, the singular values are in fact the Schmidt 
coefficients of 1^;), thus we suggest to name them generalized Schmidt coefficients. Notice 
that Yli=o ^fi^) — since the pure quantum state \z) is normalized. The number r of non- 
null generalized Schmidt coefficients is the rank of C{z). Note that the Schmidt criterion to 
decide whether a given pure quantum state is entangled or not is readily contained in this 

brmulation, since \z) is separable if and only if rc{z) = 1, which is equivalent to o'i{z) = 

27|. 

In matrix notation this all is elementarily formulated as follows: 
Proposition IV. 1. Let \z){z\ G PS'i)7Vi(g)Af2 ■ Then 

1. \z){z\ = {F{z)(^G*{z)) \a){a\ {F^ (z) ® G\z)) , where F{z) (G{z)) is the left (nght) 
matrix of the GSVD of G{z) and \a) = Y^^Zq o-i{z)\ii) . 

2. \z){z\ G SEPj^^i^i^^ if, and only if, (Ti^z) = 0. 
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A. Full-rank states with r^^ = 1 



We will find those nonnegative A such that p = (1 — A)Mi ® M2 + A|z)(2;| G SEPiVi(g)7V2, 
i.e. we will compute A*. Since Vp = Ni ■ N2, Mi and M2 are also full-rank matrices, hence 
there always exist nonsingular matrices Fi and F2 such that Mj = FiF^ 



3J]. We have the 



following partial result: 

Proposition IV.2. Let p = (1 - A)Mi (8)M2 + A|;z) (^1 with M, = FiF} . Then p E SEPn^^n^ 
if, and only if, p = (1 - A)^^ + \\cr){a\ G SEPn^(^n2, where \(r) = Xli=o 0-^(2:) |m) is built 
with the generalized Schmidt coefficients of C{z) with respect to [M^^,M2^)- 

The first key result is the following necessary and sufficient separability criterion for p: 

Proposition IV. 3. Let p = (1 — A)^^i^ + \\a){a\. Then p E S'^Ptv^^atj if, and only if 



A < X*(z) 



1 



l+NiN2ao{z)ai{z) ' 



Proof. The detailed proof is rather technical and cumbersome, thus we have moved some 
technical lemmas to appendix |B] for clarity's sake. We make use of the partial transpose, 



whose positivity is a necessary condition of separability {5,6]. By lemma [K4l p^^ is positive 
if, and only if, A < A* (2) = i+jv-^jv^Jg(^)g-^(^) • Thus we must find a criterion to decide for 
which A E [0,X*{z)] we can conclude that p is separable. Since it is clearly separable for 
A = 0, if we prove the separability for A = X*{z), then from the convexity of PSDAf-^^ATj and 
hence of SEPjv^0Ar2 our result follows. The separability of (1 — X*{z))^^^ + A*(z)|cr)(cr| is 
proved in lemma [K6l □ 

This result is interesting in itself even without relation to the whole approach. Firstly, 
the state p is a more realistic mathematical description of the experimentally prepared pure 
quantum state \a) (hence of any pure state because there always exist unitary matrices U 
and such that p = (1 - A)^^ + A|^)(^| = U^V* (^{1 - X)^-^ + X\a){a\^ f/t ® \/* for 
every since it is ineludible to introduce some noise in the state preparation procedure. 
Thus states like p are relevant per se; they are indeed Werner states in arbitrary dimensions 



351] . Secondly, this result will allow us to generalize the sufficient and necessary condition of 



separability already proved for some generalized Werner states 35 
realm (see below). 



36| in the multipartite 
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Merging both results we have a first sufficient and necessary criterion of separabihty for 
bipartite states of arbitrary dimensions of the form p = (1 — A)Mi ^ M2 + X\z){z\: 

Result 1. Let p = (1 - A)Mi O M2 + X\z){z\, with Mi = FiP} , tm, = Ni, i = 1,2, and 
A G [0, 1]. Then p is separable if, and only if, A < \*{z), where \*{z) = i^jViAf2fTo(^)o-i (z) ^'^^ 
ai{z) are the generalized Schmidt coefficients of \z) with respect to (Mj""*^, Mg"^). 



B. Full-rank states with rEp > 1 



Following along the same line as above, we focus on those cases in which the decomposition 
(jl]) yields a matrix Ep such that r^j^ > 1. We arrive at the ffist incomplete result of this 
approach. 

Result 2. Let p = (1 - A)Mi ® M2 + XEp e mt PSDn,(^sN2 where Mi = FiP} , Mi ® M2 € 
int PSDn^^n2, ^ e [0,1] and Ep = Ylk=i'=^k\ek){ek\- Then p G SEPni(8N2 if ^ < A(-Ep) = 
' ^^ere A*(efe) = TTjv^jVa^feirai) • 

Proof. Let p^ = ^k'' and p = p|a=a(£;„)- Given that Y.k=i ^fe = 1 (since tr(£'„) = 1), it is 
immediately to prove that 

K 

p = ^Pfe [(1 - A*(efc))Mi ® M2 + A*(efc)|efc)(efe|] . 

/c=l 

Now this a convex combination of separable matrices, hence it is separable. Thus if p 
is separable when A = and when A = X{Ep), then from the convexity of SEPat^^atj "we 
conclude that p is separable whenever A < X{Ep). □ 

Regretfully this is an incomplete result because the condition is not sufficient, as the 
following counterexample shows. Let p = (1 - X)^^ + A (||01)(01| + where 
|^+) = -1= (|00) + 1 11)). Then X{Ep) = i, while the PPT condition states that it IS separa- 
ble for all A < ^. 

This is better illustrated in a geometric fashion (see figured]). Result [2] does not achieve 
to detect those separable states p = (1 — X)Cp + XEp with A G (|, ^] but only those within 
the "polygonal" approximation to the cone section generated by Cp, |ei)(ei| and |e2)(e2|. 
Thus far we do not have a clear hint of how to overcome this difficulty arising from the 



10 



2 X]j=l 

e2)(e2| = |^'+)(^'+| 

FIG. 1. Bidimensional section of both the cone SEP2®2 (shaded) and the cone PSD2(g)2 iUustratmg 
result ([2]) through a counterexample for its sufficiency (not scaled). 

highly non-polyhedral character of SEPAr-^07V2- 

However we leave as an open question the following: given a full-rank state p, does 
there always exist a full-rank product matrix Cp = Mi M2 such that the rank of Ep is 
1, i.e. = 1? In the light of result [1] an affirmative answer would solve completely the 
separability problem for full-rank states. On the contrary, a negative answer would entail 
important consequences to result [D the "polygonal" approximation must be overcome. 

V. QUANTUM STATES IN PRODUCT FACES 

In this section we extend the preceding results to those quantum states not having full 
rank. Notice that this approach is independent of the dimensions of the factors in the tensor 
product structure and in particular that the decomposition (jlj) is valid for any face of 
the cone PSDat^^atj. We need the orthogonal complementation of a matrix: 

Definition V.l. Let F = YUZo (A '^'^ N -dimensional square matrix with rank r. 
We define the orthogonal complementation of F , and denote it by F , as the matrix F = 
J2i=o where the vectors {|c^)}i=r,...,Af-r spans the orthogonal subspace 

to R{F). 

We agree to define F = F whenever r = N. Notice that there always exists F~^. This 
definition allows us to write the following identity: 
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ei)(ei| = |01)(01| 




Proposition V.l. Let M = FF\ with tm = r, be an N- dimensional matrix such that 
M = ^21=0 \^k){'>TT'k\, where the kets \mk) are unnormalized. Then 

where denotes the m— dimensional square null matrix and © denotes the direct sum of 
matrices. 

Proof. Since we can write F = X]i=o \^k){k\, we liave F = J2l=o l"^fc)(^l + Sil^^ 

and ttie result follows readily. □ 

We can formulate the following 

Result 3. Let p = (1 - A)Mi (g) Ms + \Ep with Mi = FiF] , tm, = Vi, i = 1, 2, = ri ■ rg, 
A G [0, 1] and Ep = Ylk=i ^k\(^k){ek\- Let Fi he the orthogonal complementation of matrix Fi, 
i = l,2. Then p G SEPn^^n,^ if \ < \{Ep) = ^r^^ziEZ' ^^^^^ ^"(^fc) = i+T,r2aXk)'^iie^) 
o'iick) are the generalized Schmidt coefficients of |efc) with respect to (^(FiPl)^^, (Fs-FJ^)^^^. 
If K = 1 the condition is also sufficient. 

Proof. The proof runs along similar lines to that of previous sections. We can decompose 



p = {l-X)Mi(^M2+\\z){z\ = Fi®F2 

(5) 

where \a) = Y7j=o'^ji^)\jj)' ^ji.^) being the generalized Schmidt coefficients of \z) with re- 
spect to {{FiFI)-^,{F2FI)-^Y Using the same decomposition as in appendix |Bl then 

(1 - \)— + X\(j){a\ IS separable if, and only if, A < A^^ [z) = 

1 

l+rir2a-o{z)cTi(z) ' 

In those cases in which Ep = ^^Li Cfc|efc) (efc| the proof is exactly the same as that of 
result El In these cases, the condition is only necessary. □ 

This result only applies to those quantum states belonging to faces which contain at least 
a product matrix of the form Mi © M2, which are those faces Fs parametrized by tensor 
product linear subspaces, i.e. faces of the form Eg = -F^^^^j. For this reason result [3] is the 
second incomplete result of this approach. 

12 



A) 



(In 



) © (1^2 ©G^JV-ra, 



rir2 
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VI. PLAUSIBLE GENERALIZATIONS TO MULTIPARTITE SYSTEMS 



Since the preceding results are independent of the dimensions of the quantum subsystems, 
they can be extended to the multipartite realm to detect genuine entanglement in the so- 
called generalized Werner states: 

Result 4. Let p = (1 - A)Mi ® ■ ■ ■ ® M„, + X\z){z\, where Mi = FiPj and tm, = Vi, 
i = 1, . . . , n. Then p is genuinely entangled if, and only if, 

^> ^ (z) - 1 + (JY^^^ n) minfc=i,...,„_i{(To(2;^^-^*l^*+^-^")o^i(-s''-'''''+^-'") ' ^^"^ 
where cr^. (2;*^ ■■■*") are the generalized Schmidt coefficients of \z) with respect to 
(0i=i ^ (S^iLi i-^j-^i)'^) viewed as a bipartite quantum state of systems ii. . 

and the rest i^+i .. .in- 
Proof. The result readily follows after applying result [1] to each bipartition of state \z), since 



to be genuinely entangled it cannot be biseparable for any bipartition 



37|. □ 



This generalizes previous results in the literature [35|, ISGj. However, notice that we can 
only conclude about the genuine entanglement of multipartite quantum states and not about 
their separability, since in this case it is necessary to distinguish between full separability 



and m— order separability 



37|. This will be undertaken elsewhere. 



As an illustration of this result let us focus upon three-partite systems and, in particular, 
upon s^tes of the form p = (1 - A)| + A|GHZ)(GHZ|, where |GHZ) is the well-known GHZ 
. Then it is immediate to conclude that p is genuinely entangled if, and only if. 



state 



A > = |. Analogously, if p = (1 - A)| + A|W)(W|, where |W) is the W state |26|, p 



is genuinely entangled if, and only if, A > j^T' ~ — sTf • 

Furthermore we see under this approach a natural connection of the SLOCC classification 
of pure quantum states with the question of separability of mixed states. Let |\[') and \z) 
be connected through the invertible local operations (ILOs) Gj, i = l,...,n, i.e. \z) = 



Gi®---®G„|^') (see |26|). Then the multipartite state p = (1 - A)Mi O ■ ■ ■ (g) M„ + A|;z)(z|, 
where Mj = FiF^ and tm, = ^i, is genuinely entangled if, and only if, the state 

Ti . . . Tn 
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is genuinely entangled, which will be the case whenever 

1 

A > 



where with relaxed notation - denote the generalized Schmidt coefficients of 

Fi'^Gi'^- ■ ■i^F~^Gn\'^)- In this case the canonical states of each class should be chosen 
with the additional criterion of easing the computation of these generalized Schmidt coefic- 
cients. In this sense, the states with permutation symmetry are excellent candidates [s^. To 

illustrate this point consider a n-qubit W state \Wn) = [|10 . . . 0) + |01 . . . 0) H h |00 . . . 

1 12" 



The state (1 — A)-|^ + A|VFri)(W^n| will be genuinely entangled if, and only if, A > 
Another example arises when considering the generalized n-qudit GHZ state \GHZlf') = 
^ Eto N • • •^)- The state (1 - A)^ + \\GHZl^^){GHZl^ \ will be genuinely entangled 
if, and only if, A > y+J^^- These results allow us to investigate the d-dimensional mul- 
tipartite genuine entanglement as the number of qudits increases: we must focus on the 
asymptotic behaviour of g{n) = rf" minfc=i... „_i{(Tq^'"**'*''^^'"*"(t]^^"'*'°'*''+^"'*"}. For instance, if 
g{n) = 0(1), then the amount of entanglement (in the sense of the amount of permitted 
noise in the pure state preparation) does not change with n. We conjecture that indeed the 
quantity g{n) can be the basis for a multiqudit entanglement measure. 

VII. CONCLUSIONS 

We can summarize the preceding results in the following partial separability check: 
Summary 1. Let p be a mixed bipartite quantum state of dimensions Ni and N2: 
i. It can always be decomposed as 



p = (1 - \)Gp + A 



K 



^efc|efc)(efc| 



k=l 



(7) 



where Gp is an arbitrary matrix such that = Tp, {|efc)} are linear independent 
vectors and K < Vp. 

a. Let p = {l-\)Mi®M2 + \ Y^k=i'^k\ek){ek\ , where Mi = FiF- , i = 1,2 and ri = tm,- 
Let Fi be the orthogonally complemented matrix of Fi, i = 1,2. Then if 

1 

Y.k=i [1 + rir2cro(efc)cri(efc)] 



A < ^ — r— ^' (8) 
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where crj{ek) are the generalized Schmidt coefficients of \ek) with respect to 
(^2^1)"^), the quantum state p is separable. 

Hi. In the same conditions, if K = 1, the condition is also sufficient. 

Reminding the cone structure of the set of unnormahzed mixed quantum states, the 
previous results show two hmitations. On the one hand, resuhs (ii) and (iii) only apply to 
those states beloging to faces with product structure = F^-^^^s^. On the other hand, 
when the universal decomposition ([7]) is such that K > 1, the condition is only necessary. 

Despite these limitations, the approach is not exhausted, since many other choices for 
the matrix Cp exist. We believe that the present limitations arise from the simple choice 
of Cp as a tensor product matrix and that a more complex structure will overcome both of 
them. This is currently under investigation. 

The independence of this approach with respect to the dimensions allows us to give the 
first steps to extend the results to multipartite systems. And in particular, result (iii) can 
be extended to detect genuine entanglement in generalized Werner states: 

Summary 2. Let p = (1 - A)Mi ® ■ ■ ■ (g) M„ + \\z){z\, where Mi = FiF} , i = 1, . . . ,n and 
Ti = TMi- Let Fi be the orthogonally complemented matrix of Fi, i = l,...,n. Then p is 
genuinely entangled if, and only if, 

X > I (9) 

1 + (nr=i ^i) minfc=i,...,„_i{(To(2;*i-*fcl*'=+i--'")fTi(z*i-*'=l*'=+i-*") ' ^ 

where cr^- (2*^ ■■■*") are the generalized Schmidt coefficients of \z) with respect to 

^1=1-1 j=ik+l 

viewed as a bipartite quantum state of systems ii . . .i^ and the rest i^+i ■ ■ - in- 

Condition dH]) only detects genuine entanglement and does not resolve between full sepa- 
rability and m— order separability. The main reason lies on the fact that decomposition ([7]) 
is thought to detect bipartite entanglement. Thus to distinguish between the different sorts 
of multipartite entanglement, we must adapt the choices of both Cp and Ep to this situation. 
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This will be undergone elsewhere. However, notice that the approach is not exhausted once 
more. 



Finally a natural connection between the multipartite separability problem and the pure 
state classification under SLOCC arises, since the entanglement of any state admitting the 
choice Cp = Mi ® ■ ■ ■ ® M„ is equivalent to that of p = (1 — A) '^^^'"j^" + XHi (8) ■ ■ ■ ® 
■ ■ ■ ®H\, where is the canonical state of the particular SLOCC class. An 
appropiate choice of will drive us to an easier computation of the generalized Schmidt 
coefficients of ifi ® ■ ■ ■ (8> if„|\E'). 
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Appendix A: The generalized singular value decomposition 



The singular value decomposition (SVD) is a standard tool in matrix analysis (cf. e. 



3J]) by which any rectangular matrix M G A^m,n(C) can be written as 



M = UDV^ = Y,^k\uk){vkl (Al) 

k=l 

where ai > cr2 > ■ ■ • > (X-r > are the singular values, the matrix U , whose columns 
i = l,...,m, are the so-called left singular vectors, is unitary and the matrix V, whose 
columns \vj), j = l,...,n, are the so-called right singular vectors, is also unitary. The 
integer r is indeed the rank of the matrix M. 



The generalized singular value decomposition (GSVD) is a possible extension of this result 



[see e.g. 



40|) which states that, given two m— and n— dimensional positive definite matrices 



A and B, any rectangular matrix M G A^m,n(C) can be written as 
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M = UDV^ = Y,at''\ut'){vt\ (A2) 



k=l 

A,B A,B ^ ^ AB 



where a^' > ^2 ' > ■ ■ • > fx,^'^ > are the generahzed singular values with respect to 
[A, B), the matrix tj, whose columns \uf'^), i = 1, . . . ,m, are the so-called left generalized 



singular vectors, is A-unitary j^l] and the matrix V", whose columns \vf''"), j = l,...,n 



A,B\ 



are the so-called right singular vectors, is 5- unitary 42j. The integer r is again the rank of 



the matrix M. In the text we have dropped out the specification {A, B) to ease the notation. 

Notice that the difference between the SVD and the GSVD is just the scalar product 
of the rows and the columns of M under which the standard proof of the decomposition is 
estabhshed (see {m]). 

Appendix B: Lemmas 

We gather some technical lemmas necessary to prove the preceding results. 

Lemma B.l. Let oj^ , where = 0, 1, . . . , n — 1, he the nth complex roots of the unit. Then 

^ , f Vp such that ^ # N, 
I n Vp such that f G N. 



k=0 



Proof. Elementary 43|]. □ 



Lemma B.2. Let / = {0, . . . , r — 1} and lo = e^. Let uq = 1, Ui = u*, U2 = OJ^ , = w**^, 
. . . ,ujm = CJ^*-*^™,. . . , m > 1, be complex phases where the complex conjugation appears only 
in the odd indices. Then for all i,j,p, q & I we have 

{1 ifi = j^p = q or i=pj = q , 

(B2) 
ijj'^ijvq otherwise, where rriijpq G Z — {0}. 

Proof. The expression ujiUj*uj*ujq, with uji = 1, u* , u"^ , u^* , co^ , u^^* , . . . , reduces to 1 only in 
the stated cases, that is, ii i = j,p = q or i = p, j = q. In the rest of cases, it is (almost) 
evident that the expression reduces to a power of u or oi u*. □ 

Lemma B.3. Let ctq > cxi > ■ • ■ > a^-i > be such that YllZo ^1 = ^- Then we have 
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r-1 



r-1 



k=0 



k=0 



Both are elementary. 



(B3) 



Proof. This is equivalent to prove the following two numeric inequalities for alH = 0, . . . , r- 
1: 



□ 



Lemma B.4. Let Q = (1 - A)^^ + \\a){a\, where \a) = EI=o ^^N^)- Then Q^^ > z/, 
and only if, X < j^j^ 



croai 



Proof. The partial transpose of Q is given by 



Q 



r-1 



r—1 T—1 



-1 r-1 



i=0 j=0 



=0 j=0 



min(Afi-l,Af2-l) 



where \sij) = + and a^j = ^ — Now the system of vectors {|m)}™q 

lis. A Ir; •U°'^''^^l"^'^2"^h ) ik -\ I \yori=mm{iVi-l,iV2-l)+l,...,max{JVi-l,iV2-l) . ^.^.thouor 

i\Sij), \a'ij)u,j=o u IPjjA P«j//jor i=o,...,min(JVi-i,Af2-i) ortuouor- 

mal basis in C^^ C^^. Thus the identity matrix IniN2 admits a spectral decomposition 
in terms of the unidimensional orthogonal projectors onto these vectors. Summing up all 
the terms we have the spectral decomposition of Q^^ with eigenvalues of the alternative 
forms + Xcrf, -^^^ + XcTiO-j and — Xa^aj, all of them will be nonnegative provided 
— XuiUj > for all i ^ j- Since the generalized Schmidt coefficients constitutes a 
nonin creasing sequence, then A > ^ 



l+7ViAr2(Tocri ■ 



□ 



u 



For the next result we need to define M = maxijpg^j{mijpg} and Uq = M +1, where 



m. 



ijpq 



and I are defined as in lemma IB.2t r < min(A'^i — 1, Ai"2 — 1) and are normalized and in 
non- increasing order. Define u 



2-iTi 

e "0 . 
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Lemma B.5. We have the identity 



r-l 



1=0 



\ 



k=0 



r—1 r— 1 



/ .-1 \ 



J i=0 j=0 



no 



\^k){uk\^\ul){ull (B4) 



fc=0 



where \uk) = ^i^^l^'^l^) and \ul) = J2U ^f^^^^l^) ■ 
Proof. Taking into account the definitions we have 



'r-l 



r-l 



5^fT, ) \Uk){Uk\ ® \ul){ul\ = J2 ® 

ijpq=0 



J=0 



which after summing up in the powers of k yields 

r—l /no— 1 



*iP9=0 \ fc=o 



(B5) 



(B6) 



Taking into account lemmas IB.ll and IB. 21 for the choice of the base no of the roots, we 
will have 44 1 



r-l 



r-l 



ij=0 



ij=0 

Dividing by uq we finally get the desired result. 



Lemma B.6. With the same definitions as above, we have 

r r 

i=i j=i 



Proof. Firstly, since In,n2 = E^I'o ^ ® Si=o ' 1^)01 ^^'^ ^^^^^ 



(B7) 



□ 



(B8) 



r-l r—l 



i=0 j=0 
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r-1 N2-I Ni-1 r-1 Ni-1 N2-I 

2=0 i='' i=0 i=r j=r 

r-1 r-1 r-1 \ 

^0^1 5^ ® |J)(J| + ® • (B9) 

ij=0 i=0 j=0 / 

Since the three first matrices in the right hand side are clearly separable, we will focus 
on the last one, which we rewrite as 



r-1 



r-1 



r—1 r—1 



j=0 



1=0 



1=0 j=0 



In view of lemma IB. 31 we can write 



r—1 / r—1 \ r—1 r—1 



i=0 



r-1 



i=0 



r-1 



\ 



\ 



j=0 



i=0 j=0 



r—1 r—1 



i=0 j=0 



where the difference is clearly a separable matrix. Now, this is exactly the matrix in lemma 
IB.5t which is separable. 
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